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ABSTRACT 

Using published luminosity and molecular gas profiles of the late-stage mergers NGC 3921, NGC 7252 
and Arp 220, we examine the expected luminosity profiles of the evolved merger remnants, especially in 
light of the massive CO complexes that are observed in their nuclei. For NGC 3921 and NGC 7252 we 
predict that the resulting luminosity profiles will be characterized by an r 1 / 4 law. In view of previous 
optical work on these systems, it seems likely that they will evolve into normal ellipticals as regards their 
optical properties. Due to a much higher central molecular column density, Arp 220 might not evolve 
such a "seamless" light profile. We conclude that ultraluminous infrared mergers such as Arp 220 either 
evolve into ellipticals with anomalous luminosity profiles, or do not produce many low-mass stars out of 
their molecular gas complexes. 

Subject headings: galaxies: individual (NGC 3921, NGC 7252, Arp 220) — galaxies: interactions — 
galaxies: evolution — ISM: molecules — infrared: galaxies 



1. INTRODUCTION 

The merger hypothesis for elliptical galaxy fo rmation, 
as pu t forth by Toomre & Toomre (1972; see also Toomrc 
1977), posits that two spiral galaxies can fall together un- 
der their mutual gravitational attraction, eventually evolv- 
ing into an elliptical-like remnant. An early objection to 
this hypothesis was that the cores of ellipticals are too 
dense t o result from the dissipati o nlcss merging of two 
spirals (Carlberg 1986 , Gunn 1987, Hcrnquist, Spergel & 
Hcyl 1993 J. An obvious solution to this objection is to 
include a dissipative (gaseous) co mponent in the progeni- 
tors ( for other sol u tions, see e.g. Veeraraghavan & White 
1985, Barnes 1988, Lake 1989). Numerical experiments in- 
cluding such a component readily showed that gaseous dis- 
sipation can efficient ly drive large amounts of ma t erial into 
the central region s ( Ncgropontc fc White 1983| , Noguch] 
& Ishibashi 1986|, IBarnes & Hcrnquist 1991). Indeed, it 



was considered a great success for early hydrodynamical 
work to be able to reproduce dense knots of gas within 
the central regions of simulated merger remnants, similar 
to the gas concentrations observed in IR. luminous mergers 
(Barnes & Hcrnquist 1991, Sanders et al. 1988b). Since 
the inferred central gas mass densities in the observed gas 
knots are comparable to the stellar mass density seen in the 
cores of normal ellipticals (~ 10 2 M Q pc -3 ), it seems natu- 
ral that t hey could be the seed of a hig h surface brightness 
remnant ( Kormcndy fc Sanders 1992] ). 

Subsequent numerical work by Mihos & Hernquist 
(1994; hereafter MH94) finds that the dissipative response 
of the simulated gas component is so efficient that the 
resulting mass profiles of the simulated remnants are un- 
like those seen in normal ellipticals. In particular, ensu- 



ing starformation leaves behind a dense stellar core whose 
surface density profile does not join smoothly onto the 
de Vaucoleurs r 1 / 4 profile of the pre-existing stellar pop- 
ulation. Instead, the profiles exhibit a "spike" at small 
radii, with a suggested increase in surface brightness by 
factors ^100. While the predicted break in the mass den- 
sity profile occurs at spatial scales comparable to the grav- 
itational softening length of the simulations, making the 
precise slope somewhat questionable, the conclusion that 
the profiles should exhibit a clear break was considered 
firm (MH94). This prediction, if confirmed, offers a means 
to constrain the frequency of highly dissipative mergers in 
the past by searching for their fossil remnants in the cores 
of nearby ellipticals. However, the numerical formalisms 
used in MH94 to model gaseous dissipation, star forma- 
tion, and energy injection back into the ISM from mas- 
sive stars and SNe ("feedback") are necessarily ad hoc in 
nature. As such, these predictions should be viewed as 
preliminary (as Mihos & Hernquist themselves note). 

We investigate this question observationally by convert- 
ing the observed gas column densities of on-going or late- 
stage mergers into optical surface brightness by assuming a 
stellar mass-to-light ratio appropriate for an evolved pop- 
ulation. This light is added to the observed luminosity 
profile after allowing it to age passively. The resulting lu- 
minosity profile is examined for anomalous features such 
as the sharp break predicted by MH94. 

We conduct this experiment with the late-stage mergers 
NGC 3921 and NGC 7252 and with the ultraluminous in- 
frared (ULIR) merger Arp 220. These systems were chosen 
because they have been observed in the CO (1-0) molecu- 
lar line transition with resolutions (full width at half max- 
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2.5" (|Yun fc Hibbard 19994 |Wang et alj 
1997[ ). The resulting spatial radial 



resolution (300-400 pc; H = 75 km s _1 Mpc -1 ) is similar 
to the hydrodynamical smoothing length used in MH94 (~ 
350 pc) 2 , indicating that the molecular line observations 
have sufficient resolution to resolve the types of mass con- 
centrations found in the simulations. The molecular gas 
surface densities of each of these systems are plotted in 
Figure [j], converted from CO fluxes by adopting a con- 
version factor of Ag, /I cq = 3 x 10 20 cm~ 2 (Kkms -1 ) -1 
(Young & Scoville 1991). Detailed studies on each of these 
systems, which fully discuss their status as late stage merg- 
ers, c an be found in Schwcizcr (1996) and Hibbard fc van 



plot shows that the observed gas densities in NGC 3921 
and NGC 7252, although high, are not high enough to 
significantly affect the present luminosity profiles. The 
profile of NGC 3921 is basically indistinguishable from an 
r 1 / 4 profile. The profile of NGC 7252 does show a slight 
rise at small radii, but not the clear break predicted by 
MH94. Therefore the resulting luminosity profiles of these 
remnants are now and should remain fairly typical of nor- 
mal elliptical galaxies, and the conclusions of MH94 are 
not applicable to all mergers of gas-rich galaxies. Since 
both of t hese systems also obe y the Faber- Jackson rela- 
tionship ( [Lake fc Dressier 1986[ ) and NGC 7252 falls upon 
the fundamental plane defined by normal ellipticals ( [Hib 



Gork qm (1996) for NGC 3921; Schweizer (1982) and Hib- bard et al. 1994 [Hibbard 1995|) , we conclude that at least 



bard et al. (1994) for NGC 7252; and Scoville et al. (1997) 
for Arp 220. 

2. RESULTS 

2.1. NGC 3921 & NGC 7252 

For these moderately evolved merger remnants (ages of 
0.5-1 Gyr since their tidal tails were laun ched, Hibbard| 
et~aL 1994[ |Hibbard fc van Gorkom 1996] ), the observed 



gas and luminosity profiles are used to predict the expected 
luminosity profile of a 2 Gyr old remnant. We assume that 
all of the molecular gas is turned into stars at the same 
radii, adopting an exponentially declining starformation 
history. The present luminosity profile is allowed to fade 
due to passive aging effects, and the final luminosity pro- 
file is the sum of these two populations. The molecular gas 
profiles plotted in Fig. [j] are converted into gas mass den- 
sities by multiplying by a factor of 1.36 to take into con- 
sideration the expected contribution of Helium. Optical 
luminosity profiles have been obtained by Schweizer (1982, 
1996), Whitmore et al. (1993), and Hibbard et al. (1994), 
showing them to be well fitted by an r 1 / 4 profile over all 
radii, with no apparent luminosity spikes. 

The gas surface densities (T, gas in M pc~ 2 ) are con- 
verted to optical surface brightnesses (hb) by dividing 
by the stellar mass-to-light ratio (M*/Lb) expected for 
a 2 Gyr old population. We adopt the stellar mass-to- 
light ratios given by de Jong (1995, Table 1 of ch. 4), 
which were derived from the population synthesis mod- 
els of Bruzual & Chariot (1993) for an exponentially de- 
clining star formation history, a Salpeter IMF, and Solar 
metallicity (M*/L B = 0.82 MqL" 1 at 2 Gyr). Noting 
that 1 L pc~ 2 corresponds to [Ib = 27.06 mag arcsec -2 
(adopting Mb.q = +5.48), the conversion from gas surface 
density to optical surface brightness is given by hb{i") = 
27.06 mag arcsec -2 - 2.5 x log[Y. gas {r)/(M*/L B )}- 

The luminosity profiles of the evolved re mnants are 
estimated from the observed -B-band profiles ( |Schwciz"ci 
1996, Hibbard et al. 1994), allowing for a fading of +1 ma: 



arcsec z in the _B- band over the next 2 Gyr (Bruzual 



Chariot 1993, Schweizer 1996), and adding in the expected 
contribution of the population formed from the molecular 
gas, calculated as above. We emphasize that this should 
favor the production of a luminous post-merger popula- 
tion, since it assumes the that none of the molecular gas 
is lost to stellar winds or SNe and the adopted IMF favors 
the production of many long-lived low-mass stars. 

The results of this exercise are plotted in Figure ^. This 



some mergers of gas-rich systems can evolve into normal 
elliptical galaxies as far as their optical properties are con- 
cerned. 

2.2. Arp 220 

Since Arp 220 is an extremely dusty object, its opti- 
cal luminosity profile is poorly suited for a similar analy- 
sis. Instead, we use a luminosity profile measured in the 
near-infrared, where the dust obscuration is an order of 
magnitude less severe. Arp 220 was recently observed 
with camera 2 of NICMOS aboard the HST (Scoville et 
al. 1998), and we use the resulting AT— band luminosity 
profile, kindly made available by N. Scoville. Since Arp 
220 is presently undergoing a massive starburst, the fading 
factor is much less certain than for the already evolved sys- 
tems treated above, and depends sensitively on what frac- 
tion of the current light is contributed by recently formed 
stars. We adopt a situation biased towards the production 
of a discrepant luminosity profile by assuming that the en- 
tire population was pre-existing, converting the observed 
K— band profile to an evolved B— band profile by adopting 
a B — K color o f 4, appropriate for a 10 Gyr old population 
(de Jong 1995). The contribution due to the population 
formed from the molecular disk is calculated exactly as 
before. The resulting profile is shown in Figure |[ 

This figure shows that Arp 220 is predicted to evolve 
a luminosity profile with a noticeable rise at small radii. 
This is due to the peak in the molecular gas surface den- 
sity at radii less than 0.5 kpc (Fig. |lj). We conclude that 
Arp 220 has the potential to evolve a similar feature in its 
luminosity profile, if Indeed all of the current molecular gas 
is converted into stars. However, the expected rise of ~2 
mag arcsec 2 in surface brightness (a factor of ~6) is con- 
siderably lower than the two orders of magnitude increase 
predicted by the simulations (see Fig. 1 of MH94). 

3. DISCUSSION 

From the above exercise, we conclude that neither NGC 
3921 nor NGC 7252 are expected to show a significant de- 
viation in their luminosity profiles, and that the maximum 
rise expected for Arp 220 is considerably lower than the 
two orders of magnitude increase predicted by the simula- 
tions of MH94. We conclude that the numerical formalisms 
adopted in the simulations to treat the gas and star forma- 
tion are incomplete. Mihos & Hernquist enumerate var- 
ious possible shortcomings of their code. For example, 
their star-formation criterion is extrapolated from studies 
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of quiescent disk galaxies, and may not apply to violent 
starbursts. Perhaps most importantly, their simulations 
fail to reproduce t he gas outflows seen in ULIR galaxies 
"superwinds" e.g. Heckman, Armus fc Miley 199Cl|, Heck 



man, ^chncrt fc Armus 1993Q , suggesting that the numcr 
ical treatment of feedback is inadequate. 

In spite of these results, it is still interesting that under 
some conditions there might be an observational signature 
of a past merging event in the light profile of the remnant. 
The question is for which mergers might this be the c ase? 
Since Eg., is tight ly correlated with IR luminosity ( Yun| 
& Hil |>bard 1999bV we infer that only the ultraluminous 
IR galaxies retain the possibility to evolve into ellipticals 
with a central rise in their luminosity profiles. While such 
profiles arc not typical of ellipticals in general, they are 
not unheard of. For example, ~ 10% of the the Nuker 
sample profiles presented by Byun et al. (1996) show such 
anomalous cores (e.g., NGC 1331, NGC 4239). It is there- 
fore possible that such systems evolve from ultraluminous 
IR galaxies. This can be tested by careful "galactic ar- 
chaeology" in such syst ems to search for signatures of a 
past merger ev ent (e.g., 3chweizer & Seitzer 1992, Maliii| 
& Ha flley 1997j ). 

However, it is not a foregone conclusion that systems 
like Arp 220 will evolve anomalous profiles. This system 
presently hosts a ver y powerful expanding "superwind" 
(Heckman ct al. 1996), which may be able to eject a signif- 
icant fraction of the c old gas in a "mass-loaded flow" (e.g. 
Heckman ct al. 1999]). Such winds are common in ULIR 
galax ies flHcckman, Armus fc Miley 1990 , Heckman, Lchn- 
ert fc Armus 1993]). Another related possibility is that 



the IM F may be biased towards massive stars (i. e. "top 
heavy" , [Young et al. 198^, |5coville fc Soifer 199l|). Such 



mass functions will leave far fewer stellar remnants than 
the IMF adopted here. A third possibility is that the stan- 
dard Galactic CO-to-H2 conversion factor is inappropriate 
for ULIR galaxies, and that the high gas surface densities 
derived from CO observations (and thus the result ing stel- 
lar lu minos i ty profile) may be over-e stimated (see Downes 
et al. |1993| , [Bryant fc Scoville 1996[ ). 

Some support for the idea that central gas cores may 
be depleted by the starburst is given by a population syn- 
thesis model of NGC 7252, wh ich suggests that it expe- 
rienc ed an IR luminous phase (Fritze-von Alvenslebcn & 
Gerhard 1994). While the current radial distribution of 
molecular gas in NGC 7252 is flat and lacks the central 
core seen in Arp 220, it appears to connect smoothly with 
that of Arp 220 in Fig. [j]. Therefore one may speculate 
that NGC 7252 did indeed have a radial gas density pro- 
file much like Arp 220 but has since lost the high density 
gas core as a result of prodigious massive star formation 
and/or superwind blowout. However, the burst parame- 
ters are not strongly constrained by the available obser- 
vations, and a weaker burst spread over a longer period 
may also be allowed ( Fritze-von Alvenslebcn fc GcrharT] 
19940 Further insight into this question could be obtained 
by constraining the past star formation history in other 
evolved merger remnants. 

In conclusion, a comparison of the peak molecu- 



lar column densities and optical surface brightnesses in 
NGC 3921 and NGC 7252 suggests that some mergers 
between gas-rich disks will evolve into elliptical-like rem- 
nants with typical luminosity profiles, even considering 
their present central gas supply. For ULIR galaxies like 
Arp 220 the case is less clear. Such systems will either 
produce an excess of light at small radii, as seen in a small 
number of ellipticals, or require some process such as mass- 
loaded galactic winds or a top-heavy IMF to deplete the 
central gas supply without leaving too many evolved stars. 
If the latter possibility can be excluded, then the frequency 
of such profiles may be used to constrain the number of 
early type systems formed via ULIR mergers. 3 

4. SUMMARY 

• Even under assumptions that favor the production of 
a luminous post-merger population, the dense molec- 
ular gas complexes found in the centers of NGC 3291 
and NGC 7252 should not significantly alter their 
luminosity profile s, which are alr e ady typical of c l- 
liptical galaxies (Schweizcr 1982, [Bchwcizcr 199G ). 
Since these systems also obey the F aber- Jackson re- 
lationship ( |Lake fc Dressier 1986| ) and NGC 7252 
falls upon the fundamental plane defined by nor mal 
ellipticals (|Hibbard et al. 1994 [Hibbard 19~95[), it 



appears that at least some mergers of gas-rich sys- 
tems can evolve into normal elliptical galaxies as far 
as their optical properties are concerned. 

• The dense molecular gas complex found in the cen- 
ter of Arp 220 may result in a moderate rise in the 
remnants' luminosity profile at small radii. Since the 
molecular gas column density is a tig ht function of 
IR luminosity ( Yun fc Hibbard 1999b| ) , we conclude 
that this condition may apply to all of the ultra- 
luminous infrared galaxies. However, this does not 
preclude a merger origin for elliptical galaxies since 
(1) About 10% of the Nuker sample ellipticals (Byun 
et al. 1996) show such rises in their radial light pro- 
files, and (2) it is possible that much of the gas in 
such systems is blown into intergalactic space by the 
mass-loa ded superwinds found emanating from such 
objects (Heckman et al. 1999 , Heckman, Lehnert & 
Armus 19§3fr . 



• The maximum expected rise in the luminosity pro- 
files are considerably lower than the orders of magni- 
tude increase predicted by the simulations (MH94). 
We therefore suggest that the numerical formalisms 
adopted in the simulations to treat the gas and star 
formation are incomplete. 
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Fig. 1. — Radial surface molecular gas mass density profiles (in Mq pc 2 ) derived from the integrated CO flux maps are plotted as a function 
of radius. The CO fluxes have been turned into mass densities by adopting a conversion factor of Nh 2 /I C o = 3 x 10 20 cm- 2 (Kkms- 1 )- 1 . 
The data are from Yun & Hibbard (1999a) for NGC 3921 (dotted line); from Wang et al. (1991) for NGC 7252 (solid line); and from Scoville 
et al. (1997) for Arp 220 (dashed line). 
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Luminosity Profiles of Merger Remnants 
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Fig. 2. — Predicted luminosity profiles as a function of the fourth-root of the radius for the aged remnants of NGC 3921, NGC 7252, 
and Arp 220. Dotted lines: expected contribution to the B— band profile by the present stellar component, after aging by 2 Gyr. For NGC 
3921 and NGC 7252 this is calculated from the presently observed B— band luminosity profiles (taken from the data in Schweizer 1996 and 
Hibbard et al. 1994) and fading by 1 mag arcsec - 2 . For Arp 220 we use the observed K— band luminosity profile (from Scoville et al. 1998) 
adopting B — K = 4 for a 10 Gyr old population. Dashed lines: expected contribution to the B— band profile from the expected post-merger 
population if all the molecular gas is turned into stars. Solid lines: the expected B— band luminosity profiles after adding in the aged stellar 
and post-merger populations. No significant deviation from a pure r 1 / 4 law expected for NGC 3921 or NGC 7252. Arp 220 does not show 
such a seamless profile, exhibiting a slight rise in at small radii. 



